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DETAILED COMPUTATIONAL PROCEDURE FOR DESIGN OF CASCADE BLADES WITH PRE- 
SCRIBED VELOCITY DISTRIBUTIONS IN COMPRESSIBLE POTENTIAL FLOWS : 


By Сковає В. CosTELLO, Ковевт L. Cumuinas, and Јонх T. SiNNETTE, Jr. 


SUMMARY 


А detailed step-by-step computational outline is presented for 
the design of two-dimensional cascade blades haring a prescribed 
relocity distribution on the blade in a potential flow of the usual 
compressible fluid. This outline is based on the assumption 
that the magnitude of the velocity in the flow of the usual com- 
pressible nonviscous fluid is proportional to the magnitude of 
the relocity in the flow of a compressible nonriscous fluid with 
Linear pressure-volume relation. The computational procedure 
includes sereral ways of adjusting the prescribed velocity to sat- 
isfy restrictions imposed by the method. Tables of coefficients 
are gicen for evaluating the necessary integrals, including the 
determination of the harmonie conjugate function. Numerical 
examples are included. 


INTRODUCTION 


A method for computing blade profiles with prescribed 
velocity distributions based on the assumption that the pres- 
sure-volume relation is linear is presented in reference 1. 
The method uses the prescribed velocity distribution and 
compatible free-stream conditions to determine a mapping 
function that transforms an incompressible flow about the 
unit circle into an exact compressible flow, with linear pres- 
sure-volume relation, about a cascade of blades having the 
desired velocity distribution. 

In order to apply this method to the design of a cascade 
with a given velocity distribution end free-stream conditions 
in a flow of an “actual” fluid (the usual compressible non- 
viscous fluid), the relation between the actual fluid and the 
fluid with the linear pressure-volume relation must be approx- 
imated so thet the required velocity distribution and free- 
stream conditions for the second fluid may be determined. 
In this investigation, which was made at the NACA Lewis 
laboratory, the magnitudes of the velocities in the two fluids 
are assumed to be proportional and the constant of propor- 
tionality is determined by the continuity equation using the 
same upstream and downstream flow angles for the two fluids. 

By use of this relation, a computational procedure was 
developed to obtain the blade profile with & minimum of 
effort. "The procedure, presented herein, includes the adjust- 
ment of the prescribed velocity distribution to satisfy the 
restrictions on the mapping function and the numerical 
computation of the harmonic conjugate function. 


A, B, Cs, ... 


SYMBOLS 
The following symbols are used in this report: 
constants 
constants 


function of 8 defined by equation (18) 
spacing of cascade blades 
function of 0 


final values of Z(e**) 

Imaginary part 

distorted velocity defined by equation (22) 

modified .K (8) 

constant determined by range of potential 

Mach number 

constant determined by trailing-edge angle of 
blade 

dimensionless velocity on blade in direction of 
increasing & 

magnitude of dimensionless velocity (ratio of 
velocity to stagnation velocity of sound) for 
flow with linear pressure-volume relation 

final value of q after modification of ReH (et) 

real part 

change in ReH (6%) 

ratio of velocities (determined by equation (4)) 

аге length along blade measured counterclock- 
wise from tail (figs. 1 and 2) 

are length along final blade shape 

magnitude of dimensionless velocity (ratio of 
velocity to stagnation velocity of sound) 
for actual fluid 

velocity on unit circle (incompressible flow) 

complex variable in caseade plane 

angle of velocity in compressible flow (meas- 
ured from positive x-axis) 


circulation (positive counterclockwise) 
ratio of specific heats 


difference between Re H(e*) and parabola at 
half-interval point 


including trailing-edge angle of blade 


circle angle (incompressible flow plane) pos- 
itive counterclockwise 


1 Supersedes NACA TN 2281, “Detailed Computational Procedure for Design of Cascade Blades With Prescribed Velocity Distributions in Compressible Potential Flows" by George R. 


Costello, Robert L. Cummings, and John T. Sinnette, Jr., 1951. 
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À auxiliary function defined by equation (15) 
р density 

T variable of integration 
T velocity potentia] 
Subscripts: 

с compressible flow 

1 incompressible flow 

n leading edge (nose) 

і trailing edge (tail) 

1 upstream 

2 downstream 


RELATION BETWEEN FLOWS 


In the present problem of designing a cascade, the pre- 
scribed conditions are the upstream velocity we’, the 
downstream velocity u,e%2, and the magnitude of the di- 
mensionless velocity on the blade as a function of the аге 
length u=«u(s) (see figs. 1 and 2). The upstream and down- 
stream velocities are related bv the continuity equation 


1 


us) Uz, сово;  (l) 


ele 


= -1 | = 
(1-255 79} и, сов a-(1 -— 1 


Consequently, only three of the quantities 25, М», ол, and 
аз muy be assigned, and the fourth is determined by equa- 
tion (1). 

‚ In utilizing the method of reference 1 to design the cas- 
саде, the prescribed conditions are employed to determine 
the upstream and downstream velocities and the velocity 
distribution on the blade for the fluid with the linear 
pressure-volume relation by assuming that the magnitudes 
of the velocities for the two flows are proportional; that is, 
$: q=ru (2) 
The constant of proportionality is determined from the 
additional assumption that the flows have the same free- 
stream directions by substitution of equation (2) in the 
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FIGURE 1.—Caseade fn z-plane. 
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FIGURE 2.—Presoribed velocity distribution on easeude blade, 


continuity relation for the flow with the linear pressure- 
volume relation: 


| | А | 
-——— 0 COS ,=—==——— d, COS о, (3) 
Vr ds vita’ 

, Ue? COS? a,—"u;? COS? ay ^ 
. Au us (cos a ceos? aj) | 


© (4) 


—- 


Hence, the upstream velocity ge’, the downstream ve- 
locity фе‘, and the velocity distribution on the blade 


9=94(8) are now determined for the flow. with the linear 
' pressure-volume relation. 


Other approximations to the relation between the two 
flows could be used but, with the present cascade inter- 
pretation, the upstream and downstream flow angles and the 
spacing are the same for the two fows. The two densily- 
velocity relations are also in good agreement, as shown in 
figure 3 where the density-velocity relations are plotted for 
u,—0.700 and u:=0.500. When the upstream and down- 
stream velocities are equal, the density-velocity relation, as 


given by this interpretation, is the same as the density- 


velocity relation in the Kármán-Tsien approximation. 
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FiGURE 3.—Comparison of density-velocity relations for иін 0. 700 and 14-0. 500. 
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NUMERICAL PROCEDURE 


After the free-stream conditions and velocity distribution 
on the blade are determined for the fow with the linear 
pressure-volume relation, the method of reference 1 is used 
to obtain the blade shape. In this method, the given 
eonditions are used to select an incompressible flow about 
the unit circle and to determine the mapping function that 
transforms this flow into a compressible flow with the circle 
transforming into the desired blade. 
for determining the circle flow and the mapping are given in 
the following paragraphs. 


FLOW IN CIRCLE PLAXE 


The velocity potential ¢.(s) and the circulation Г, on the 
blade are obtained by integrating the prescribed velocity 
distribution g(s). For this purpose, it is convenient to 
define Q(s) by 
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Епис ве 4.— Velocity and potential distribution for fluid with linear pressure-volume relation. 


A detailed procedure. 


where the total are length has been taken equal to 27. 
These integrals are evaluated by Simpson's rule using suf- 
ficient points so as to obtain accuraté values of ¢,(s). The 
values of ф.(з) should be plotted on the graph with Q(s) 
(fig. 4) in order to obtain q easily as a function of the circle 
angle later in the computation. Fhe spaemg d of the 
cascade is given by 


Г, 
Фон (7) 
qı SIN а — G2 SIN c | 
' From equation (35) of reference 1, the meompressible 
potential function for points on the unit eirele may be 


written 


tan 6 


cos 8 
TY 


cosh Ё 








ei(0) = —2 Re Аташ + Tm A+iJm B)tan 


sin @ 


(Im A—Im В) ап! cen Oy d (8) 





where D is chosen to make the potential zero at the tail 
stagnation point б,, the angle convention is—;«tan ahr 


< and ап! c Is taken in the same quadrant and the 


same direction as 0. "The values of He 4, Im A, Re B, and 
Im B are determined from the free-stream conditions and 
spacing by 





Re а= паза (9) 

Im Á——41-4 gj Re Atana (10) 
Re B=—Re А (11) 
ImB=y1+ ¢.' Re Atanas (12) 


The value of the constant & is determined by the condition 
that the potential range on the cirele is equal to the potential 
range on the blade: 

e«(2) T 8) = v8 +2х) E pila) (1 3) 
The proper value of k is computed by assuming а value for & 


and computing the stagnation points 8, and 6,, which are 
the roots of the equation 


—(m AtIm В sinh k 


sin (8-Е)= 2 Re Asecr (14) 
where А is given by 
A—tan^! um d га алаг (15) 
T T 
E 
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and then determining the value of ¢,(@,+27)—¢,(6,). This 
process is repeated for different assumed values of Ё to ob- 
tain a plot of ¢,(6,+27)—¢,(6,) as а function of k. Thevalueof 
k satisfying equation (18) is then obtained by interpolation. 
A computational outline of this process for use with a desk 
calculator is given in table I, in which the usual column op- 
erations are written as lines. 

When k has been determined, the flow about the circle 
is known with the value of the potential on the circle given 
by equation (8) and the velocity on the circle given by 


4 Re А cosh k sec XT. | P 
"(o оре зр | sin (83- X) —sin 4+») | (16) 
A computational outline for v(0) and фч(@) based on 80 points 
is given in table IT, lines 1 to 19. 


DETERMINATION OF MAPPING FUNCTION 


The mapping function (equation (62), reference 1) may be 
written in the form | m | 
4 g;?* —v(8)* (2 cosh 2k —2 cos 29) Anio +з tani Er de 


4 g, €*4/2 cosh 2k —2 eos 20 


z= 


(17) 
where 


g= [2 —2 cos (6,— 8)] niie [Re Cfe#)+(2 cosh 2k—2 cos 20) Re H(ett)] 


£2=n n (H ") т етно cosh 2k — 2 cos 2@)Im H (ef?) 


and 


€ 
n=] —- 
т 


where є 18 the trailing-edge angle of the blade in radians. 
The mapping function is completely determined when 
Re Cte"), Im C(e"), Re H(e*), and Im H(e**) are known. 

Determination of Re C(e®) and Im C(e**).—The function 
Се?) is defined by equation (44) (reference 1) and may be 
written in the form 


d Aes з gi) | 


(еі нод, ARG) 


1 mr É— kei B 
gie O peg, 


Ce) = 5(1 ecu) In 
| T 8) 


Hence Re C(e*) and Im C(e*) are obtained by taking the real 
and imaginary parts of equation (18): 


ReC(e'*)—5| (C+ C)--&(0,— C) cos веб; C.) sin J 
(19) 


ImC(e)=4 | (C+ С) еке, С) cos à— e (C,— C) sin J 
(20) 


where 
C;=I1n iby rae HIS RAFT т? Ае [е ®(2 cosh k— 
1 


2 cos кт 


2 соз M 


Im A =j —sin 0, 
Re n ш (= cos 8, 


тър me 8 
=: ича =! ma- E " г 


The quadrants for Ше are tangent terms in C; and C, are se- 
lected by considering the numerator and the denominator of 
the argument as signed quantities and choosing a quadrant 
in which the sine has the same sign as the numerator and 
the cosine has the same sign as the denominator. For con- 
venience, these angles are taken as positive and the resulting 
values of C; and Cj are changed by multiples of 2x until Cs 
and |C,| are each less than 2r. 

Determination of Re Н (e'5).—In order to obtain Re H (e*), 
the prescribed velocity along the blade must. be determined 
as a function of the circle angle g=q(@). This relation is 
obtained from the equality of the potentials at. corresponding 
points in the two planes. Because Q(s) and ¢,(s) are plotted 
on the same graph (fig. 4), the magnitude of the prescribed 
velocity at а given circle angle 0 may be found directly by 
reading the absolute value of Q (g—lQ[) at the abscissa 
where e,(8) equals the calculated value of ¢;(6) (table IT, 
line 19). The prescribed velocity is thus obtained as a 
function of the circle angle g=q(@). 

From equation (48) (reference 1) 


Pee ГУУ 
L, К(9) [2 — 2 cos (0, — 6)]"/* Ве C(e uid 
| 2 cosh 2k — 2 cos 20 





C,—tan^! 


In 
Re Н (е) 


(21 ) 
where 
. .2g(0 . 


yr 000008 


Lines 20 to 36 of table II show the detailed computation for 

obtaining Re H(e*), including the evaluation of Re C(e*). 
Adjustment of Re H(e").—Restrictions on the mapping 

function require that Re Н (е) satisfy the conditions 


К(0)= 


f ” ВеН(еч) de=0 (23) 
f Re H(e")sin 9 40--0 (24) 
p Re Н (ef?) cos 8 d8—0 (25) 
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These integrals are evaluated numerically (in the manner 
developed by Dr. Glenn Н. Peebles, while at Douglas Air- 
craft Co., Inc., in his work on isolated airfoils in incom- 
pressible flows) by using a parabolic variation of Re H(e") 
between three consecutive points and integrating analyti- 
cally. The point spacing used in the examples of this 
report is an equal spacing of 80 points (designated whole 
points and numbered 0, 1, 2, . . . 79) around the circle 
beginning with point 0 at é=—180° and continuing around 
the circle at 4.5? intervals. Where Re H(e*) has large 
fluctuations (usually in the neighborhood of the stagnation 
points), the midpoints of these intervals (called half-points) 
are also used. The coefficients for evaluating the integral in 
equation (23) are the Simpson's rule coefficients. The 
coefficients for the integrals in equations (24) and (25), 
which are obtained by integrating analytically the product 
of the parabolas and sin @ or cos 6, are given in tables ПІ to 
VI. When half-interval spacing is used, it should begin and 
end at even-numbered points of the original 80 points and 
the integration coefficient to be used at the beginning (or 
ending) even-numbered point is one-half the sum of the 
whole-point coefficient and the half-point coefficient at that 
point. 

If the values of Re H(e*) from equation (21) do not satisfy 
equations (23), (24), and (25), the velocity distribution 
chosen is incompatible with the chosen free-stream conditions 
and must be modified to make these integrals zero. The 
most desirable method of modification will depend on the 
magnitude of the integrals in equations (23), (24), and (25) 
and on the features of the original prescribed conditions that 
ure to be preserved. Usually the original free-stream con- 
ditions are retained and the velocity distribution is modified. 
When the integrals are large, however, changing the circle 
flow may be desirable in order to minimize the changes in 
the essential characteristics of the original velocity distribu- 
tion (such as limits on maximum velocity, diffusion rate, 
апа so forth). 


Modification of the circle flow provides essentially one 
additional degree of freedom, as can readily be seen from 
equations (9), (10), (11), and (12). The strength of the 
singularities А and B are determined by the free-stream con- 
ditions and the spacing. The range of potentiel is deter- 
mined by the location of the singularities, that is, by the 
constant k. Only the ratio of tbe potential range to the 
spacing is important, however, because changing both to- 
gether merely changes the dimensions of the cascade by a 
scale factor. This additional degree of freedom may be con- 
veniently represented by the ratio of are lengths of the upper 
and lower surfaces. Consequently, this ratio has a large in- 
fluence on the size of the integrals in equations (23), (24), 
and (25). The selection of the ratio may be based on the 
ratio for a blade having a similar velocity distribution; or, 
iu some cases, it may be advantageous to try several ratios 
of arc lengths and roughly approximate the integrals using 
only 20 points and from these results select the proper ratio 
to minimize these integrals. The final adjustment to reduce 
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these integrals to zero may then be obtained with & fixed 
circle flow by adjustment of Re Н (е'). 

When adjustments are being made in Re H(e*), the change 
in profile are length and profile velocity 4(0) produced by 
these adjustments must be considered. The local arc length 
ds corresponding to d6 will be decreased in the same ratio 
that the velocity is increased in order to maintain the same 
potential. The change in arc length is automatically ac- 
counted for in the final integration for the blade coordinates. 
The change in the velocity produced by changing Re H(e*) 
by an amount A(Re Н) is indicated by the change in the 
"distorted" velocity A(0), which is related to g({@} by equa- 
tion (22) (fig. 5). 
given by 


K (8) — K (8) е 2 cosh 25-2 cos 20 A (Re Н) (26) 


Making a change in Re H(e*) therefore has the effect of mul- 
tiplying the distorted velocity by а factor. This factor varies 
widely with the eircle angle and produces the smallest change 
in K(é@) at 0— +r ог 0—0 and the largest change at вб-- + x/2. 
Hence, for a given percentage change in K(8), the greatest 
values of A(Re Н) occur near 0= +r or 6=0. The effect of 
k on these values of A(Re Н) for а given ratio of КК) 
is shown in figure 6. 

Consequently, for small values of the mtegrals in equations 
(23), (24), and (25), complete adjustment can frequently be 
made by merely changing the slope of the distorted velocity 
K(@) at the stagnation points. The changes in the integrals 
produced by this change in slope are given by 


г A(Re Ну de 


|" A(Re H) sin 0 48 


IL A(Re H) cos 0 d0 
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FIGURE §.—Distorted velocity КВ) as function of (0). 
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where A(ReH) is determined from equation (26). These 
changes in slope are effective in eliminating or reducing 
S Re Н(е“) 29 and f Re H(e'*) cos 6 do. (бее figs. 6 and 7.) 

Another method for adjusting Re H(e*) is to use а multi- 
plier on the entire velocity distribution; that is, 


K(0) —e-^K() 
then 
Ate) e nate, fly wos ote dads 
2 cosh 2k —2 cos 28 
and 


Га (Re H)d8— 5. 


f A(Re Н) sin 0 d@=0 


Г A(Re Н) cos 0 d6—0 


Changes in the sine and cosine integrals may be obtained 
by using different multipliers over the upper and lower 
surfaces: 


K()—K(9 е 
É (9) — K(0) вто» 


8, 0X6, 


0, X 0 8,--2x 
then 


А (Re Н)- 8, «8X0, 


2 cosh 2k—2 cos 20 | '— — 


ыз ДЫ... e 
A (Re H) 2 cosh 2k —2 cos 26 ва S 0 S Or 28 


and 2 
т = ay _ tan 0 Vs 
Га Re н)4в=(„ sinh 2Ё ven tanh А Је, m 
tan 6 важ 


-1 
n tanh k/a, 


(тїшїр 
2 sinh 2k 





Ы . SEN Gi -1 cos 8 
| ден) віп 0 d@= (= соат tanh сову 3 Aa 


Qs _, cos Ê в, аж 
(5 cosh k балів" E), 


x = G1 —1 sin 0 Зп 
(7 ден) cos 0 40 (iar SEE tan sink): + 
з Sin 0 \ fet 


an: sinh Ё /e, 


(zs 
4 sinb k . 





These integrals may be easily evaluated because the inverse 
trigonometric functions all enter into the computation of 
¢:(@) and have been calculated in the determination of Й. 
A possibility for making small adjustments in Re Н (е) is 
to use a multiplier that is a simple function of 6, such as 
K= 


КІ е «9 oos a+c аїр б) 


then 
_ a+b cos @+c sin 8 
A (fie H)—2 cosh 2k —2 cos 20 

and 


L med) do aie 2k 


f” A (ReH sin eg tr ib D... 


т аза 07 (1 —tanh E) | 
Г ден) cos 0 de= зу р pa 
Various combinations of these methods for adjustment of 
Re H(e*) have been used in the illustrative examples of this 
report. 

After the adjustments in Re Н (е) have been made, the 
final values should be checked in equations (23), (24), and 
(25), as indicated in table П, lines 42, 43, and 44. At this 
time, Ё (8) should be computed in order to determine whether 
É(8)« 2, as required by part (d) of equation (5), reference 1. 
Н É(0) does not satisfy this inequality, a different modifi- 


cation of Re H(e*) is necessary. 
ШЖ ЕИ ШЕШ 
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FIGURE 6.—A(Re 7T) as function ої k and в for ratio of distorted velocitles ДФ / А (б) 0.050. 
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FIGURE 7.—Products A(Re H) oos 0 and A(Re H) sin 0 for k&0.20 and for E($/ (Фф) 0.050. 
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The final values of Re H(e*) апа ¢(@) should be plotted in 
order to insure that no excessive fluctuations have been 
introduced. At this time, q(8) should be compared with the 
original q(0) to make certain thet the velocity has not been 
changed beyond acceptable limits. Because the potential 
range is fixed, these changes will change the arc length. Н 
desired, 2(@) may be obtained as a function of the new arc 
length 5 by 


a= q[&(8)] = @(8) 


~~ f" ре (6)| 
ra= |, AG de 


and 





0, <8 <4,42x 


The values of : E аге given in table И, line 50. 





Computation of Im H(e'*).—When Re Неї?) satisfying all 
requirements has been obtained, the computation of the 
conjugate function Im Не? is done by direct numerical 
integration of Poisson's integral 





Im He =5- f Вее") cot $=" dr 

on the SO basie points using an extension of the method 
developed by Dr. Glenn H. Peebles. The integration is 
accomplished by replacing Re (е9) by 40 parabolas on the 
80 basie points and integrating analytically the product of 
the parabolas and the cotangent term. The two sets of 
80 coefficients so obtained for the integration are given in 
tables VII and VIII. (See appendix A for derivation.) 
The set in table VII is used for the computation of the con- 
Jugate at the end points of the parabolas; that is, at the even- 
numbered points 0, 2, 4, . . . 78. The set in table VIII 
is used for the computation of the conjugate at the odd- 
numbered points 1, 3,5, . . . 79. Thus, to obtain 7 m H(e* 
at one of the basic points, for example, point 7, the value of 
Re Неї at this point is multiplied by the first coefficient 
(0.000000) in table ҮШ, the next value of Re Ef (e?) (at 
point 8) is multiplied by the next coefficient (0.412368), 
and so forth, and the sum of these 80 products is the desired 
value of Im Й (е4) at that point. Hence to obtain Im (е) 
at the 80 points requires 80 such accumulative multiplica- 
tions. This computation is done very efficiently on an elec- 
tronic calculating punch using only 415 cards and taking 
approximately 2 hours, including the time for key punching 
and verifying. 

When Re He) has large fluctuations so that some of the 
80 basic intervals are divided by half-points, the preceding 
values of Іт Й (єї? must be corrected to take into account 
the difference between the values of Re (ef) and the basic 
parabolas at these half-points. The correction to be added 
to the calculation is the harmonic conjugate of these differ- 
ences. The coefficients for computing the value of the con- 
jugate (at the 80 basic points) of these differences are given 
in table IX. (See appendix В.) This computation has been 
arranged on two concentric disks—the coefficients are carried 
by the outer disk and the differences (denoted by 5) are 
entered on the inner disk in the proper places (fig. 8). It is 

212483—54——-8 





FIGURE 8.—W beel for computing correction 2 Im Hes). Indicated setting gives correction 
at po 15. 


easy to see where these differences will make a significant 
contribution to Zm Ё (еі) and to obtain the contribution by 
accumulative multiplication. This correction is then added 
to the calculated results. 


COMPUTATION OF BLADE COORDINATES 
After the corrected values of the conjugate have been 
obtained, the blade coordinates are given by equation (17), 
which for convenienee of computation may be written in the 
parametric form 
sin 28 


m d BOI т =f - 
а 3 (9) cos СЕ (0) —tan m ME 8 (21) 








2 | . т in 20 
y= (^ ZO an| otita 0 tan? 275 де em 
where 
(n (HN Mn C (e*) -- (2 cosh 2 — 


2 cos 28) Im E (еч) 


These integrals are evaluated by Simpson’s rule. In order 
to obtain sufficient accuracy, use of at least the same number 
of points as were used in the evaluation of Re H(e*) is 
advised. Because the values for Im Й (e'*) were calculated 
at the 80 basic points only, interpolation of g»(8) is necessary 
when half-points are used. The values of r(@) and y(@) are 
given in table II, lines 70 and 71. 


EXAMPLES 


Several examples have been computed to illustrate some 
of the variety of conditions to which the method may be 
readily applied. No attempt has been made to pick the 
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best velocity distribution or to obtain the most desirable blade 
shape, as to do so would lead to many considerations far 
beyond the scope of the present investigation. 

In these examples, the adjustment of the prescribed veloc- 
ity altered the arc length slightly. For uniformity, figures 
9 to 12, which show the final velocity distribution and blade 
shape, were therefore scaled to give an arc length of 27. 

Example 1.—For this example of a cascade blade with 
low turning and a cusped trailing edge ae the free-stream 
conditions were taken аз .. .  ..— "m MP" 


440.488 oe,—170? | 


140.478 205—180? .. 


The velocity distribution prescribed on the blade was the 
velocity distribution for an isolated Joukowski airfoil. Ad- 
justments to Re Н (е) altered this distribution somewhat 
and the resulting blade (which is very similar to a Joukowski 
airfoil) and velocity distribution are shown in figure 9. 


Example 2.—For an exemple of ап impulse-type blade. 


with rounded leading and trailing edges (e—3), the free- 
stream velocities were assigned the values 


Uy; =. 555 
Ug=0. 555 


a, == 135° 
047225? 
and the prescribed velocity distribution on the blade was 


constant over most of the upper and lower surfaces with 
different values on the two surfaces and varied linearly 
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FIGURE $.—Final velocity sip e and blade anane for Example 1. Mı=0.500; «=0 
cusped trailing edge). 


a"i 70? 
| м,” 0468 


through the stagnation points. In this example, the expres- 
sion for r (equation (4)) is indeterminant and the value of г 
was obtained by taking the limit as а, approaches 225°. 
In this example, complete adjustment of Re H(e*) was made 
by multiplying the velocities on the upper and lower surfaces 
by constants. The resulting blade shape and velocity dis- 
tribution are shown in figure 10. A smaller nose or tail 
radius of curvature could be obtained by increasing the slope 
of the velocity distribution through the nose or the tail, 
respectively. 

Example 3.—The free-stream velocities chosen for this 
example were representative of a compressor stage with 


Vu4—0. 583 
43-0. 417 


a, — 135? 


аз = 1 55° 


In order to keep the changes in the prescribed velocity to з 
minimum, three ratios of lower surface length to upper surface 
length of 0.95, 0.90, and 0.85 were used with the same pre- 
scribed velocity distribution (but, of course, different circle 
flows) and the integrals in equations (23), (24), and (25) were 
quickly approximated using only 20 points. By use of these 
results, & ratio that would give the smallest integrals was 
chosen (& ratio of 0.93) and the blade shape with rounded 
trailing edge (e—«) was computed in the usual manner. 
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FIGURE 10. oe velocity distribution and blade shape for Example 2. —.Mi»0.573; «вх 
(rounded trailing edge). 
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The blade shape and the velocity distribution are shown in 


figure 11 (a). By use of the same circle flow and a change 


in @(@), a blade having а 10° trailing-edge angle ( ‘=T5 and 


similar velocity distribution was obtained (fig. 11 (b)). 

Example 4.—In this example of a highly loaded cascade 
blade with a rounded trailing edge (e— x), the free-stream 
velocities were 


Hy =. 579 
м„==(). 369 


a; = 135° 

a= 180° 

The prescribed velocity distribution, which was similar to 
the distribution on a Griffith airfoil, had an abrupt decrease 
in value on the upper surface for use with suction. The 
blade and velocity distribution are shown in figure 12 (э). 
By use of the same circle flow and a change im ¢(@), a blade 


having a cusped tail (є==0) and essentially the same velocity 
distribution was obtained (fig. 12 (b)). 


DISCUSSION | 


Specification of the trailing-edge angle of the blade requires 
that (в) have a zero of order ос ats=Qands=2r. Ifthe 


prescribed velocity does not go to zero in the proper manner, 


у 
а-928 
u,” 04/7 T 0, = 155° 
= ця 2583 
с 7/35* 
= х 
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S 
(а) ewx (rounded trailing edga). 


Re H(e*) will beinfiniteat6—8, Fairing Re H (e) smoothly 
through 0-8, and keeping the values finite will, however, 
insure obtainmg the specified angle at the trailing edge of the 
blade. The shape of the blade in the immediate vicinity of 
this point will depend on the velocity prescribed in this 
vicinity. If desired, the blade shape at the tail can easily 
be modified by changing Jm H(e*), which changes the angle 
of the tangent to the blade (the angle in equations (27) and 
(28)). The corresponding change in Re H(e'*) is computed 
from Poisson's integral with the constant term zero, and 
the modification in the velocity is obtained from this change 
in Re H(e*). The change in Im H(e*), denoted by A (Im Н), 
should be chosen to satisfy the following conditions: 


[` sam mae=0 
f ” A(Im Нусозё 48=0 
LH AU m H)sin 0 d@=0 
For & prescribed velocity distribution requiring only mod- 
erate changes, the blade shape can be obtained in approx- 
imately 50 computing hours using 80 points and 4-decimal 


accuracy. The time depends, however, on the degree of 
familiarity with the method and on the extent of the permis- 
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(b) e= уд (199 trailing-edge angle. 


FIGURE 11.—Final velocity distribution and blade shape for Example3. Mie0.904. 
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(a) етт (rounded trailing edge). 
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. (b) ew 0 (cusped trailing edge). 


FIGURE 12—Final velocity distribution and blade shape for Example 4. Mq: 0.599. 


sible modification to the velocity. When several examples 
have the same free-stream conditions and cascade spacing, 
the time is considerably reduced because the circle flow nee 
be eomputed only once. | 


With the calculation based on a spacing of 80 points with 
half-points around the nose and the tail, the method gives 
accurate results in all cases in which the parameter & is not 
less than 0.10. Experience has shown that for cascades of 
moderate stagger and turning k will usually be greater than 
0.10 when the solidity is less than 1. In examples 1, 2, 3, 


and 4, the values of Ё were 0.2600, 0.1100, 0.2008, and 0.2851 
and the solidities were 0.91, 1.01, 0.92, and 0.74, respectively. 
In applying the method to other cases, a finer point spacing, 
which would require new coefficients for the integrations, 
should be used. 


Lewis FLIGHT PROPULSION LABORATORY 
NATIONAL ÁpvISORY COMMITTEE FOR ÁERONAUTICS 
CLEVELAND, Оно, August 28, 1950 
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APPENDIX A 
COEFFICIENTS FOR DETERMINING CONJUGATE FUNCTION 


The coefficients for determining Im H(e*) are obtained 
by considering the contributions of one of the parabolas to 
the conjugate at one of the basic points. The contribution 
I of the (E-F-1)th parabola through the points Re Б (еә), 
Re He*x+1), and Re Hí(e'z42) to the conjugate at any 
point 0— 6, is given by 





Bo 
Tax | [alana -00— Ores) + 
Re Н (вічні) cot TE = de 
(Al) 
where 


_ He Hitet8e+2}—2 Re Н(е'їзьч1)-Ь Де H(eiftzi) 
о}? 


Re Н(е'%%+т)у— Re H (e*t) 


= э} 


and А із the length of the basic interval. By substituting 
bayi On +lh 


in equation (А1} and expanding in series, equation (А1) can 
be integrated to give 


1= В, Re Не 6-7» х) S, Re H(e! (Ac HS) З 
T, Re H (ei о tl) 


2) res) asso (2) ш ні 


- ) Б ОРІ 


where (Рог [0, 1) 
ІГ ЄКМ, 1 
ве з (5) Tis 
о (y+ cha nL : 
пути - б 





. +1 
sin EN h 
= In aT SE ИН ТАНЕ e ) 
sin 5 h 
А 2N 4 17 
T=} — TE )+тв(5) +l) T2535 2) +-- |+ 
2 (VDP РИГА VL з ува (2+ 


туга ру .. yf 


and 





fe 2) Е (tan 5) JS 
LOGRO PS 3 
(tan (m5) | г (ноз) JH 








(27+5) 5 — tarn EE M 
г (заз) 5) Е 
ей (09 — 


| (2) | 1 1 (ња 5) 
DV Leryn вна 3 У 


1 44(tan 3) 


23(tan 5) 








(2j 4- 5) 15 (2j-7) 105 
Ач 
! 563(tan 2) | 
(2749) 1575 
When /=0, | 
Уся 0 
1 1/ÀAN 1 hẹ 2 FRY 
=R=] 1-5 (5) 03) -a 2) - 
aama) 
42,525 2) 
When /—1, 
s, У ДЕ M Юю ^ MA LL LL e) 
+ «M2 36 1350 26,460 425,250 


1/h*,9h*, А 4h§ 
i= 18 078 $820 212,628 6287 7" a) 
and В; is undefined because the СЕ becomes infinite. 
The final integration coefficients, denoted by No, №, 
Nz, . . . are obtained from the preceding values of £j, S; 
and Т... When 6, is an end point of a parabola (m is even), 
then 


М,=0 
` N: =S; 
зна (A9) 
Na, 1 Sg 
Na Tuc Hai 


where №, is determined by Cauchy” s principal EE When 


0, is а midpoint of a parabola (m is odd), then 
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N o=S)=0 The coefficients for 80 points given in tables VII and VIII 

Мү=Т,-+ В; were obtained from equations (A2) and (АЗ) using А==2+х/80. 

NS The coefficients for other point spacings can be obtained from 

e (АЗ) | equations (А2) and (АЗ) by using the proper value of А. 

n Because of the symmetry of the coefficients, only half of 

Мис Тара Ra them need be calculated. 
Матея 
APPENDIX В 


COEFFICIENTS FOR CORRECTION OF CONJUGATE 


The coefficients used in the wheel correction are obtained 
by determining the contribution of one of the half-point 
differences to the conjugate at one of the basic points. 
Denoting the contribution of the difference. &+4 (the dif- 
ference at 0=0:+4) to the conjugate аё 06—6, by J, then 


1 fi 45 8-- On 
I=. Е (£—5. 4). | eot 5 d 


which, on substituting 





h 
0,44 = On ta tlh 


becomes 


р рені T 
н [- | е-( nta ptt) + to 


—@„ 


ЖОП teot o rim do |. (BD 





Integration of equation (B1) by series expansions gives 
(for 1250) 


{| sm 0+1) А 


L= snt Ја — y —2Un6 BGP + EQ + ice .) 
вн 
where 
an l 
4 
tan (214-1) т - ue 
and 


E tan i) | i. 1 2( tang) i) | Зо) 
ЕРА 5 / |і) (273) 





ЕТ 7145 . 
44 44 (ав). ЭЙ 563 (tan 1) | 
ОГ ^ 105 | tern» ^ 1575 
When [—0, 
I—$ S (3 35$ 33200 338680") 
— "nthe \4 288 43,200 3,386,880 


Consequently, the coefficients for the wheel correction 
denoted by АЗ, Ni+a, Мачі, . . ., аге 


_8 ( -aks узро 33803907 " 9 
М я 4 288 43,00 3,386,880 


шй рі з 
Мы т In — ур 2G + E46G5-- E3G5-- ...) 
sin = I 
2 


The coefficients in table IX, which are for 80 points, were 
obtained from equation (B2) using А = 2x/80. By use of 
other values of À in equation (B2), the coefficients for dif- 
ferent point spacings can be obtained. Because of the sym- 
metry of the coefficients, only half of them need be computed. 
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TABLE I—OUTLINE FOR COMPUTATION OF k 








Line Operation Remarks 

1 | Assume a value of k 
2 | sinh qi 
8 | cosh (1 
r WIPE 

т - in 
b Re A x (4) 
6 + tan-t (5 (0) Is ^; express In degrees; ~90°<A< 90° 
7 | sec (6) | 
8 (Im A+Im В),,(2) 

є 2Re А (7) 
р | sin- (5) 9) їв Atin; - 909 «4-8, 3909 
10 | (9) —(6 10) із ĝa 
ii 7. (10) --2(0) — 180° 11) is 6, 
22 | (11) о 
18 | sin по 
14 | соз (10 
15 | tan (10) 
16 | sin "2 
17 | cas (12 
18 | tan (12) 
19 | (16) + (2) 
20 | (13) -(2) 
21 | 617) + (3) 
22 | (14 + (3 
23 8) - (4d 
24 15) -- (4 
25 | tan-! (19) Express in radians; —3< (25) «5 
26 | tan-! (20) Express in radlans; -;«00« 
2; | tanh- (21 
28 | tanh-! 
29 | Чап! (23) Express ш radians; same quudrant and 
nas ёг r 
30 | tan-! (24) Express 2 rodians; same quadrant and 
8 85 da 


еве 


Must, equal gc(2x})—e.(t,) for correct 
value of і 


с» 
D 


Gm Aa Im om —(28)] 
Gn E (30)] 





бетер озыны 


5 


Operation 


Assign values of б at 4.5° intervals Ни 


inning at ё=— 
2d £a x 


cosh 2k— (3) 
4 Re A cosh k sec X) + (4) 


(Im .4— Im B)Xu8)— 
2Re AX(7) + 
(Im 4% m B)Xt18)+D 


(27) X (23) 
dj 


i 

КЫ RA 

Cy— Coe (10) 

C3— COe* X (11) 

СУ C)+63)+-Gij+2 
— 65) 3-kE-- (29) 


Геза 
-F 


Г (36) sin Ф dé 


- 


х 
Ѓ (38) cos # dé. 


Remarks 


ints may be added Ister to these 
80 points where necessary 


sec Х Is line 7, table I 
X is line 6, table I 


sin Is line 8 ,table 
m ep. e 8 а 


- 9 


Express In radíans;—7 « (16) «7 
Express ín radians; same quadrant and 


as 6 
(19) Is (A and 
sin б: 


D=—(Im A—Im В) (салі ВЕЕР 


LL 
- E LL 
2Re A tanh БЕ 


tan ё; 
(Im A-- Im B) tan: he 


Obtain by matching potentials 





(23) Is the distorted velocity ХХ?) 
Omit these lines when 2=0 


v- (2 


(23) = (23) when n=6 


Is Re Clet) 
із Re Het) 


Evaluate by Simpson's rule 
ета я by coefficients in tsbles III and 


кте by coefficients in tables IV and 


TABLE III—WHOLE-POINT COEFFICIENTS FOR SINE 


«> сю yoart СӘ N e © 


—.00821 | —.0: —. ШЗ 
--. 00620 | —. 04241 | —. 05178 
— O3 | —.08923 | ~. 10178 
—.01620 | —. HG: —. 04056 
—. 09669 | —.00609 
"Те 04936 —„ 4 
"Те 10:76 —. 
—.03081 | —. 05178 | —. (4241 
тт | —.10433 | —. 00 





INTEGRAL 





HERI 


4 





TABLE V—HALF-POINT COEFFICIENTS FOR SINE 


INTEGRAL 
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TABLE II—OUTLINE FOR COMPUTATION OF BLADE COORDINATES 


Line Operation 
40 | A(Re H) 
41 (36)-+ (10) 
42 T (41) d$ 
i 
43 * (41) sin & de 
Г 
44 f * (41) cos 8 dé 
=r 
45 — (10) X (29) 
48 e OS 
47 (18) X Q3) 
48 4—(47}2 
49 £(47) = (48) 
50 1С931-- (49) 
51 C- Сех (11) 
2 | (ce сажі 
55 [(Crt- Druid 7 (62)}+2 
54 Uneorrected Im Heit 
55 Correction for Im Ре) 
з | HG) 
вто | 100-+-180°]х5 
(29) XC 
58 | 57.29518 ((68)-+(68)} ` 
ба (ED 
бі ee 
x 61) + es 
et tan (63) 
85 | (DX(00-rF90^—(69) 
65 | sin 
Si | (50) (68) 
со x | 
-0 f ! (o8) 48 
-T 
t 
т1 (69) d$ 
TABLE iV— 


INTE 


a2 се | са спн Gera 


TABLE VI—HALF-POINT COEFFICIENTS 


WHOLE-POINT a TS 


Remarks 


Changes in Re H(eif) to make integrals 
zero; see text 


(41) is Re Heit} 
Mnst be zero 


Must be zero 


Must be zero 


(47) is Ew: must hare Eo 
(equation (5d), reference I) 
(19) is я (6) 


ії 
(50) is 


(53) Is 1m Се) 


Obtained from machine calealstlon using 
coefficients In tables VIII and IX 


Obtained from differences between Re filein) 
and раво using coefficients from 
table 


(56) is Im Hei) 
(57) =0 when a=0 


(50) Is g1(8) In degrees 


ti + 

ыа ыр Ты ЙЕ рого 

have same as (62) and (62), 
respectively 


со) Is r(8); evaluate Integral by Simpson's 


кр y (8); evaluate integral by Simpson's 
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- — -— ee 


FOR COSINE 
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TABLE VII—COEFFICIENTS FOR CONJUGATE AT EVEN- TABLE VIII —COEFFICIENTS FOR CONJUGATE AT ODD- 
NUMBERED POINTS NUMBERED POINTS 


lo. оооооо |o. 020051 Lo HR | onc 
. 636402 | . 030217 | . 0154 


«з бо =ї oh Ob» F2 с | 





мус =1 Са см нь OR 523A CO 
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TABLE РЕ БОН CORRECTION TO CON- 











| TRE = -i сард DM E XL: = mue tmi oe Ao. 
10 20 30 40 70 

0 0.010061 |o. 008014 јо. 003282 |—0. 000164 —Q. 021304 
1 ‚017188 | .007408 | 002890 | —. 000492 . 02404 
2 015596 | . : er —. 027496 
3 014222 | 1006308 | 002175 | —. 001153 —. 081974 
4 018021 | .005809 | . O01 —. 001 —. 038045 
5 011960 | .005335 | . 001488 | —.001829 —.0467 
8 011011 | . 004885 | . 001153 | —. 002175 —. 060498 
7 010156 | . 004455 | .0 a —. 085301 
8 009379 | .004042 |. —. 002890 —. 144611 
9 . 003645 | „000164 | —.003282 —. 638565 





